Disorder and the effective Mn-Mn exchange interaction in 
Ga^jMn^As diluted magnetic semiconductors 
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We perform a theoretical study, using ab initio total energy density-functional calculations, of 
the effects of disorder on the Mn — Mn exchange interactions for Gai- x Mn x As diluted magnetic 
semiconductors. For a 128 atoms supercell, we consider a variety of configurations with 2, 3 and 4 
Mn atoms, which correspond to concentrations of 3.1%, 4.7%, and 6.3%, respectively. In this way, 
the disorder is intrinsically considered in the calculations. Using a Heisenberg Hamiltonian to map 
the magnetic excitations, and ab initio total energy calculations, we obtain the effective J^ n ~ Mn , 
from first (n = 1) all the way up to sixth (n = 6) neighbors. Calculated results show a clear 
dependence in the magnitudes of the jf n_Mn with the Mn concentration x. Also, configurational 
disorder and/or clustering effects lead to large dispersions in the Mn-Mn exchange interactions, in 
the case of fixed Mn concentration. Moreover, theoretical results for the ground-state total energies 
for several configurations indicate the importance of a proper consideration of disorder in treating 
temperature and annealing effects. 

PACS numbers: 71.55.Eq, 75.30.Hx, 75.50.Pp 



I. INTRODUCTION 



The exciting possibilities of manipulating both the spin 
and the charge of the carriers in semiconductors, in such a 
way that new devices may be designed, have brought a lot 
of attention to the study of diluted magnetic semiconduc- 
tors (DMS) in the past ten years or so. Even though the 
DMS have been known for a long time?i it was the discov- 
ery of ferromagnetism in p-type (In,Mn) As systems^ that 
spurred the research in this field. This was even more so 
after the successful growth of ferromagnetic (Ga,Mn)As 
alloys. 3 This latter system has become almost a paradigm 
in the field of DMS materials. It has long been known 
that isolated MnQ a substitutional impurities give rise to 
acceptor states around 0.1 eV above the top of the va- 
lence band. Thus, the Mn atoms have a double func- 
tionality in the Gai-^Mn^As alloys: they provide both 
(i) the magnetic moments, and (ii) holes to intermediate 
the interaction between them. This somewhat simplistic 
view is much more complex than it seems at first sight. 
Ferromagnetism in Gai-^Mn^As only occurs for large 
Mn concentrations of a few percent. As a consequence, 
the acceptor levels form a band which, due to the rather 
localized character of the defect state, has a dispersion 
which is far from what would result from a free quasi- 
particle picture. Moreover, the intrinsic disorder coupled 
to this somewhat narrow band indicates that any theo- 
retical description based on an effective mass description 
should be viewed with caution. To further complicate the 
issue, in order to obtain the necessary high Mn concentra- 
tions the growth temperatures cannot be too high, which 
causes a lot of defects to be present in the samples, like 
Mn interstitials (Mn/) and arsenic antisites (Asc a )- As a 
result, the critical temperature and hole concentration, as 



a function of Mn composition, are crucially dependent on 
the details of growth conditions, 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 

In view of all these facts, it would be important to 
have a way to estimate the Mn-Mn exchange interac- 
tions (i) with as few assumptions as possible, (ii) which 
would treat the host and the Mn impurities at the same 
level of accuracy, and (iii) which would furthermore in- 
clude the effects of disorder. This approach of implicitly 
tracing out the holes degrees of freedom has been im- 
plemented in a variety of ways based on self-consistent 
methods. Van Schilfgaarde and Mryaso\*i 7 » have per- 
formed calculations of total energies, within the atomic 
spheres approximation, to extract exchange couplings, 
J's, for specific (i.e., not randomly chosen) clusters of 
closely spaced Mn ions; their results suggest a tendency 
of a decrease in \J\ when more Mn atoms are added to 
nearby sites. More recently, Xu et aL 18 used muffin-tin 
orbitals to investigate the dependence of the exchange 
coupling with the Mn-Mn distance at much larger (8.3%) 
concentrations of Mn atoms; they found a considerable 
scatter in the values of the exchange couplings. In a se- 
ries of theoretical studies, in which the effect of random- 
ness/disorder is described by the coherent-potential ap- 
proximation (CPA), Kudrnovsky et a/ii2*22 and Bcrgqvist 
et alm^ have used a tight-binding linear muffin-tin or- 
bital method, together with the magnetic force theorem, 
to study the dependence of the Mn-Mn exchange cou- 
plings and critical temperatures with the concentration 
of Mn impurities in III-V and group IV DMS. Also, Sato 
et alM*^ have used muffin-tin type potentials together 
with a KKR-CPA approach to study Curie tempera- 
tures and exchange interactions in III-V DMS. Moreover, 
Sandratskii and BrunoS 4 *^ 5 *^ have used the augmented- 
spherical-wave method within the local-density approxi- 
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mation to investigate exchange interactions, Curie tem- 
peratures and the influence of the clustering of Mn im- 
purities in (Ga,Mn)As. One should notice that the use of 
non-full potential muffin-tin-style approaches is not ad- 
equate to treat the electronic structure of covalent semi- 
conductor systems such as (Ga,Mn)As DMS. Further- 
more, we will show that disorder plays an important role 
which may not be adequately treated by simple effective- 
medium approaches such as the virtual crystal approxi- 
mation (VCA) or CPA. 

In this work we perform large supercell total energy 
calculations, based on ab initio density functional theory 
(DFT) methods. Within this approach, we treat disorder 
configurations in which the Mn atoms randomly replace 
Ga atoms. By considering two, three, and four Mn atoms 
in a supercell with 128 atoms, we cover three Mn con- 
centrations, 3.1%, 4.7%, and 6.3%, and present results 
for the effective exchange interactions, J^ n_Mn , between 
two Mn atoms which are n-th neighbors in the Ga sub- 
lattice, with 1 < n < 6. Also, in a few cases the Mn 
atoms are placed in predetermined positions, in order to 
compare the exchange coupling of two nearest-neighbor 
Mn atoms in the presence of other Mn atoms, placed at 
various separations. The present results indicate a clear 
decrease in the magnitudes of the J^f n ~ Mn with the Mn 
concentration x; from now on, the Mn-Mn superscript in 
jMn-Mn w yj om jtted, in order to simplify the nota- 
tion. 

The work is organized as follows. In Section II we 
provide a brief description of the calculational procedure 
used in the present study. Results and discussion are left 
for Section III, and Section IV summarizes our findings 
and conclusions. 



II. CALCULATIONAL METHOD 

We have performed total energy calculations based 
on the density-functional theory (DFT) within the 
generalized-gradient approximation (GGA) for the 
exchange-correlation potential, with the electron-ion in- 
teractions described using ultrasoft pseudopotentials. 27 
A plane wave expansion up to 230 eV as implemented in 
the VASP codeSS was used, together with a 128-atom fee 
supercell and 4 L-points for the Brillouin zone sampling; 
these -L-points are non-equivalent, due to the presence of 
Mn impurities. The positions of all host GaAs atoms as 
well as substitutional Mn in the supercell were relaxed 
until all the forces components were smaller than 0.02 
eV/A; our GGA lattice parameter for undoped GaAs 
turned out to be 5.74 A, which is in accordance with 
other estimates, e.g., that of Ref. |2tJ For a 128 atoms 
supercell, we consider a variety of configurations with 
2, 3 and 4 Mn atoms, corresponding to concentrations of 
3.1%, 4.7%, and 6.3%, respectively. Since calculations for 
all possible disorder configurations with more than 2 Mn 
atoms per cell is prohibitively costly in terms of computer 
time, we have considered typical configurations, as gener- 



ated through the Special Quasi-random Structures (SQS) 
algorithm^ A configuration a is generated by placing 
the Mn atoms at Ga sublattice sites (64 possible sites). 
We then calculate the pair correlation functions, up to 
the sixth- neighbor, given by: 

n m ( < r) = ^-^A m (i,i)^. (1) 

Here H m is the mth-neighbor pair correlation function, 
Z m is the number of rath-order neighbors to a site, 
A m (i, j) is 1 if sites i and j are mth-order neighbors, and 
zero otherwise; and Si is a variable taking values 0, if site 
i is occupied by Ga, and 1 if it is occupied by Mn. For a 
perfectly random (R) distribution of Mn atoms, the pair 
correlation function does not depend on m, II m (i?) = x 2 , 
where x is the Mn concentration. For a given configura- 
tion we calculate the deviation from randomness as 

6H(a) = X)( n mW - n m (i?)) 2 . (2) 

n I 

The above quantity indicates how random the a con- 
figuration is. We perform an exhaustive search over all 
possible configurations and choose to work with the ones 
with lowest 511. 

For each chosen disorder configuration, we adopt the 
following strategy within our DFT-GGA calculations. As 
an initial guess, we take all valence electrons of each Mn 
atom aligned with each other, corresponding to S = 5/2 
as expected in a d 5 configuration, and calculate the total 
energies for this configuration, as well as for an increas- 
ing number of flipped Mn total spins. The energy dif- 
ferences with respect to the aligned states, {AS}, are 
then described by an effective Heisenberg model with 
appropriate first-, second-, and so forth, up to sixth- 
nearest- neighbor Mn-Mn couplings, J n , n = 1 — 6; see 
Sec. IIIII and Appendix for details. This procedure has 
been applied^i*^ to the case of two Mn atoms in a su- 
percell with 128 sites, and we were able to infer the de- 
pendence of the effective couplings with both the Mn-Mn 
distance and direction. We have found^i^ that the cal- 
culated J n exchange couplings lead to a Mn ferromag- 
netic state, with the holes forming a relatively disper- 
sionless impurity band, and therefore that a conventional 
free-electron-like RKKY interaction should be ruled out 
as the origin of the Mn-Mn ferromagnetic coupling. 

One should notice that for two Mn atoms, if the 
spins are treated quantum mechanically the above men- 
tioned energy difference corresponds to that between the 
state with total spin 5/2 and the singlet one, leading 
to J ? 1 Q) = AE n /15. If the spins are treated classically, 

Jn — 2AE n /25 and the two approaches are entirely 
equivalent, apart from an overall multiplicative factor of 
1.2. For more than two Mn atoms, we consider classical 
spins and note that this approximation, though not cap- 
turing full details of the excitation spectra, is still able 
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to provide overall trends of the low energy magnetic ex- 
citations for a finite number of spins. 

As a final methodological comment, we note that we 
have checked for spin-orbit effects (in the case of two Mn 
atoms) through the projector augmented-wave (PAW) 
method^, and found a change from 0.29 eV to 0.24 eV 
in the total energy difference between the excited anti- 
ferromagnetic and ground state ferromagnetic Mn-spin 
alignements. Although a systematic study in this sense 
would certainly be important, this is beyond the scope of 
the present study and we have chosen to ignore spin-orbit 
effects in the total energy calculations presented in this 
work. We believe this approximation would not alter the 
general conclusions of the present study. Moreover, other 
possibilities, such as a non-collinear ferromagnetismfSiS 
have not been considered at this stage. 

III. RESULTS AND DISCUSSION 
A. Two Mn atoms 

Let us first discuss the case of two Mn substitutional 
atoms in the 128-site supercell. We considered all con- 
figurations corresponding to all inequivalent positions 
within the supercell, i.e., Mn-Mn distances varying from 
4.06 A up to If .48 A. Our total energy results yield a Mn- 
Mn ferromagnetic ground state in all cases. The relevant 
Heisenberg Hamiltonian in this case is 

H = J,iSi ■ Si+ n , (3) 

for each relative position in the supercell, where, for the 
sake of comparison with the cases of three and four Mn 
atoms (see below), S is taken as a classical spin of magni- 
tude 5/2; n is a vector connecting nth-nearest-neighbor 
Mn atoms replacing Ga atoms. The estimates for J n thus 
obtained are displayed in the second column of Table [I] 
As previously noted^i*^ the resulting Mn-Mn ferro- 
magnetic effective coupling in Gai-^Mn^As is essentially 
intermediated by the antiferromagnetic coupling of each 
Mn spin to the quasi-localized holes. Also, the observed 
non-monotonic behavior of J n should be attributed to the 
anisotropic character of the effective interaction. More- 
over, \J n \ essentially decreases^LiS with Mn-Mn separa- 
tion and vanishes above ~ 11.5 A. 



B. Three Mn atoms 

In the case of three Mn atoms in a supercell with 128 
sites, we have performed calculations for 10 different dis- 
order configurations. Figurc^shows two SQS illustrative 
configurations: in (a) the 3 Mn atoms are somewhat clus- 
tered together, whereas in (b) two are nearest neighbors 
and the third is farther apart. 

For each disorder configuration, the relevant Heisen- 
berg Hamiltonian must contemplate the possibility of in- 
teractions occurring not only amongst spins within the 



Table I: Estimates for the effective exchange coupling, J n , in 
meV, between nth-nearest-neighbor Mn spins, Si and Sj, for 
different Mn concentrations, x. In the case of 2 Mn atoms 
(x — 3.1%), J n is unique for a given n. For 3 Mn atoms 
(x — 4.7%) in a supercell with 128 sites, we have performed 
calculations for 10 different SQS disorder configurations, and 
J„ is given by the average over the configurations in which two 
Mn sites are nth-neighbors; the number of such configurations 
are shown in square brackets, and the error bars are calculated 
as standard deviation of averages. In the case of 4 Mn atoms 
(x — 6.3%), we show the results for two configurations (see 
text); note that sometimes a specific configuration would not 
accommodate the pertinent J„. 



n 


x = 3.1% 


x = 4.7% 


x = 6.3% 


x = 6.3% 


1 


-23.2 


-18.2 ± 1.5 [7] 


-12.6 


-13.0 


2 


-10.4 


-3.8 ± 1.8 [4] 




-4.7 


3 


-13.6 


-6.6 ± 2.7 [7] 


-2.8 


-6.0 


4 


-5.6 


-3.6 ± 0.8 [4] 


-4.8 




5 


-2.6 


+0.4 ± 0.7 [5] 


+0.1 


-1.3 


6 


-4.4 


-1.9 ±0.7 [2] 







supercell, but between one spin in the supercell and 
the different images in neighboring supercells (periodic 
boundary conditions effects, PBCE's). In actual fact, 
depending on the disorder configuration, the same pair 
of spins may be j-th nearest neighbors within the su- 
percell and k-th nearest neighbors when the images are 
considered. One can therefore write the Hamiltonian as 

H =} J },w i j(n)J n Sj • Sj, (4) 

n i<j 

where the my (n) are geometrical weights taking into ac- 
count PBCE's. For a given configuration, one expects 
most of the id's to vanish; also, we set w = if the 
distance between the Mn atoms is larger than 11.5 A, 
as previously established^ In the Appendix we discuss 
the Hamiltonian for the two SQS configurations of Fig. 
n We then calculate the total energies for different Mn 
spin configurations: with all spins aligned, with only one 
spin reversed, cither in site f , 2, or 3, and so forth, in- 
creasing the number of spin flips, until one has the same 
number of unknowns ( J n ) as equations (namely, the cor- 
responding energy differences with respect to the aligned 
state). 

It is instructive to lay out J n as a function of n for 
the 10 SQS realizations of disorder (three Mn atoms), as 
shown in Fig. [21 for comparison, we show the results for 
x = 3.1% in the same figure. One can see that the over- 
all trend of J n with n, observed in the case of two Mn 
spins, is maintained in this case, with the non-monotonic 
behavior still being due to effects of directionality, i.e., 
the exchange coupling depends not only on the distance 
between the pair of Mn atoms, but also on their relative 
direction with respect to the bonds of the host GaAs. 
Here we should mention that our 128-atom supercell to- 
tal energy results for the ferro- and antiferromagnetic 
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Figure 1: (Color online) A pictorial view of two possible real- 
izations of disorder for three Mn atoms in a 128-site supercell 
(x = 4.7%). Ga sites are represented by the smaller spheres, 
As sites by the middle-sized ones, and Mn atoms by the largest 
ones. For clarity, supercells are repeated along the different 
cartesian directions. The three nonequivalent Mn atoms are 
shown as different shades of gray (blue, red, and yellow in the 
color version). 

states are in overall agreement with the corresponding 
64-atom supercell total energy results of Mahadevan et 
al.wiL The corresponding average values of J n , for each n, 
are shown in the third column of Table [I] It is interest- 
ing to note that all J n decrease (in absolute value) as the 
concentration of Mn atoms increases from 3.1% to 4.7%. 
While at first sight this may seem an unusual behavior, 
one should have in mind that the effective Mn-Mn inter- 
action is hole-mediated, thus sensitive to the hole density. 

In order to assess the effects of clustering in a system- 
atic way, we have also considered non-SQS configurations 
in which two Mn atoms are first neighbors, and a third 
Mn atom is placed in positions corresponding to fifth- 
, third-, and first-neighbor of the pair: we found that 
J\ = —20.8 meV, —17.3 meV, and —8.1 meV, respec- 
tively; the extreme values are shown in Fig. [21 as filled 




1 2 3 4 5 6 
n 

Figure 2: (Color online) The nth-nearest-neighbor exchange 
coupling as a function of n for x = 3.1% and x = 4.7% 
(with data displayed for 10 SQS configurations, see text). Full 
curves are guides to the eye (for x = 4.7% the full line goes 
through average values of Jn). Filled squares for Ji corre- 
spond to extreme values obtained for the non-SQS configu- 
rations; see text. Also shown is the multiplicity of each n-th 
neighbor pair in a given direction < hkl >. 



Table II: Total energies from ferromagnetic SQS configura- 
tions, labelled from £ = 1 to 10, with respect to the total 
energy of the configuration corresponding to three nearest- 
neighbor Mn atoms clustered together. The effective Heisen- 
berg Hamiltonian can be written in the form H = aJiSi ■ 
S2 + /3JjSi ■ S3 + 7JfeS2 • S3, such that the entries in the third 
column are {i a j^k" 1 } 



l-Xh state 


E t (eV) 


{i a fk~<} 


1 


0.054 


{l^H 2 } 

{lW} 

{l^ 1 } 
{l^V} 
{lW} 

{l 1 2 1 3 1 } 
{l^ 2 } 
{2 1 4 2 6 4 } 
{3 1 4 2 5 2 } 
{2 1 3 1 5 2 } 


2 


0.077 


3 


0.085 


4 


0.090 


5 


0.091 


6 


0.103 


7 


0.125 


8 


0.158 


9 


0.176 


10 


0.227 



squares. Thus, clustering tends to weaken the magni- 
tude of the nearest- neighbor coupling. One may attribute 
this behavior as most likely resulting from the Coulomb 
repulsion between the holes, which leads to their delocal- 
ization as the Mn atoms group together, being therefore 
detrimental of their role as mediators of ferromagnetism. 

If, on the one hand, clustering tends to decrease the 
magnitude of the nearest-neighbor exchange, on the other 
hand it leads to the energetically most stable configura- 
tion; this is in agreement with recent results from calcu- 
lations restricted to pairs of transition metals^i In Ta- 
ble[n]we display the energies of calculated ferromagnetic 
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Figure 3: (Color online) Isosurfaces for the net local magneti- 
zation m(r) (see text for definition) in the case of three Mn&i 
defects [for the configurations depicted in Fig. IT^ a)]. with (a) 
all spins aligned and (b)-(d) only one flipped spin. The green 
surface corresponds to a value of + 0.005 e/A 3 , and the blue 
surface to — 0.005 e/A 3 , with e being the electron charge. 
The black (red) spheres denote the Ga (As) atoms. 



SQS configurations relative to the clustered one in which 
the three Mn atoms are first-nearest neighbors. We note 
that the SQS configurations labelled from 8 to 10, which 
have the highest total energies of the set, correspond to 
cases in which there are no first-neighbor pairs of Mn 
atoms. Since Gao.97Mno.03As is only stable at growth 
temperatures in the range 200-300 Cfi£& the scale of en- 
ergies shown in Table [H] indicates that not many config- 



urations can be thermally activated. Clearly, there are 
several other mechanisms at play - such as mobility of 
Mn atoms, possibility of trapping on interstitials, and so 
forth -, which are not included in the present approach, 
and will determine the final distribution of Mn atoms. 

Figures EJa)-(d) show the net magnetization m(r) = 
Pl(r) — Px(r), where p a is the total charge density in the 
cr-polarized channel, for three Mn atoms with all spins 
aligned and for only one flipped spin, for the configu- 
ration depicted in Fig. ^ ( a )- Note that the densities 
on the upper right and upper left corners in each figure 
are related to a Mn atom and its image in a neighboring 
supercell. Similarly to the m(r) of one^i and two Mn 
impurities^ in a supercell, near each Mn atom the lo- 
cal magnetization has a d a — like character, whereas close 
to the As neighbors, the character changes to p s — like, 
where a — (f or [) and a = (J, or |). Also, m(r) has 
a rather localized character. The flipping of spins intro- 
duce nodes on the m(r) and subtle changes mostly on the 
orientation of the p— like lobes. Close to the Mn atoms, 
however, the local magnetization is not very sensitive to 
the flips. 



C. Four Mn atoms 

For four Mn atoms, we have considered only two dis- 
order configurations, chosen according to the SQS algo- 
rithm. A Hamiltonian similar to Eq. may be written, 
with the addition of terms involving the fourth spin, hav- 
ing in mind that the interactions with spins on image sites 
are more frequent in this case. 

For instance, in one of the calculated SQS configura- 
tions, the effective Hamiltonian becomes 

H = 2J5S1 ■ S2 + J1S1 • S3 + J3S1 • S4 

2J 4 S 2 • S 3 + J 3 S 2 • S 4 + J1S3 ■ S 4 , (5) 

where the absence of a J2 second-neighbor interaction 
should be noticed. Calculations of total energies for all 
Mn spins parallel, and for the four possible single flips, 
lead to four excitation energies, from which the J„'s 
(n 2) may be inferred. Analogous considerations apply 
to the other SQS configuration. The results are shown in 
columns 4 and 5 of TableQ] One sees that the overall ten- 
dency of J n is to decrease in magnitude as n is increased, 
in a pattern similar to that for smaller concentrations, 
though the dispersion cannot be properly assessed due 
to the scarcity of data. We also note that, as in the 
case of three Mn atoms, calculations with a non-SQS 
configuration with the four Mn atoms clustered together 
indicate that clustering decreases the magnitude of the 
first-neighbor J\ exchange coupling: J\ — —6.5 meV in 
this case, which should be compared with the —12.6 meV 
and —13.0 meV values of Table [I] 
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responding to four Mn atoms clustered together as first 
nearest- neighbors . 

From Fig. 01 we see that, in most cases, the magnitudes 
of the exchange couplings decrease as the concentration 
of Mn atoms is increased. Further, this decrease may be 
quite significant; for instance, the magnitude of the aver- 
age Ji decreases by the order of 50% when one roughly 
doubles the concentration from 3.1%. We also see that 
for the configurations in which the Mn atoms are clus- 
tered together, |Ji| also decreases as x is increased. This 
overall decrease with x can be taken as numerical evi- 
dence that a steady increase in the concentration of Mn 
atoms is not sufficient to rise the critical temperature, 
since the exchange couplings will eventually be weakened. 
Clearly other effects may be playing important roles. For 
instance, within our present approach, the hole density 
is assumed to be the same as that of Mn atoms, which, 
as mentioned in the Introduction is not really the case. 
The presence of Mn interstitials and Mn-As complexes 
also need to be taken into account in order to reach a 
quantitative agreement. Nonetheless, one expects that 
the trends unveiled here are indicative of the actual ex- 
perimental situation. 

It is important to have in mind that several theoret- 
ical works have previously examined the dependence of 
the exchange couplings with the Mn-Mn separation or 
with the Mn concentrationi 18 i 19 i 20 i 21 i 22 i 23 i 24 i 25 i 26 Some 
predict an alternating sign for the exchange coupling, 
but these predictions should be taken with extreme care, 
since these theoretical calculations are based on non-full- 
potential muffin-tin-type potentials which are not reli- 
able to treat the electronic structure of covalent semi- 
conductor systems such as (Ga,Mn)As DMS. Also, dis- 
order quite certainly is not adequately taken into account 
within simple effective-medium approaches such as VCA 
or CPA, as fluctuations in the Mn positions essentially 
lead to variations in the Mn-Mn exchange-coupling pa- 
rameters, as apparent from Fig. 0] 



Figure 4: (Color online) Dependence of: (a) J\, (b) J2, (c) 
J3, and (d) J4 with the concentration of Mn atoms. For 
x = 3.1%, J n is unique for a given n. Values for the SQS 
configurations are shown as empty circles, while the filled di- 
amonds correspond to the extreme values obtained for the 
non-SQS configurations; see text. Dotted curves are guides 
to the eye through the average values of J„. 



D. The dependence of J n with the concentration 

The data in Table [Q can also be used to discuss the 
dependence of J n with x, for a given n. In Fig. 0] we 
plot Ji, J2, J3 and J4 as functions of x. For the case 
of Ji, we also show (as filled symbols) three values ob- 
tained for the non-SQS configurations: two as mentioned 
before, in the case of three Mn atoms, and the one cor- 



IV. CONCLUSIONS 

We have performed ab initio total energy density- 
functional calculations for two, three, and four substi- 
tutional Mn atoms in a 128 atoms supercell, correspond- 
ing to concentrations of 3.1%, 4.7%, and 6.3%, respec- 
tively. In this way, we have treated the host and the 
Mn impurities on equal footing. The effects of disor- 
der have been assessed at different levels of approxima- 
tions, depending on the concentration of Mn atoms: for 
x = 3.1%, all possible non-equivalent positions of the 
Mn atoms have been considered; for x = 4.7%, ten non- 
equivalent configurations have been generated through 
the SQS algorithm, while three specific ones have also 
been considered in order to discuss the effects of clus- 
tering; and, for x — 6.3%, two SQS and one non-SQS 
configurations have been investigated. While the rela- 
tion between the densities of holes and of Mn atoms is 
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one of the yet unsolved issues in the context of DMS, here 
we have assumed that each Mn atom provides one hole; 
since our results relate to general trends, they may be 
carried over to the actual experimental situation of only 
a fraction of Mn atoms contributing with holes. It is 
also interesting to note that the cut-off of 11.5 A (which 
would correspond to x ~ 0.042) imposed on the range 
of Mn-Mn exchange couplings would appear to be in di- 
rect contradiction with experimental data by Edmonds et 
al£, according to which ferromagnetism is seen for dop- 
ings as low as ~ 0.015 (where one would have essentially 
no compensation) . Since the site percolation threshold 40 
for FCC lattices is 0.20, for the Ga FCC sublattice in 
(Ga,Mn)As, the concentration cut-off for ferromagnetic 
order would be of the order of 0.20 x 0.042 = 0.0084, i.e., 
x ~ 0.84%, indicating that there is no contradiction with 
the measurements of Edmonds et al£ 

We have focused mainly on the effective exchange in- 
teraction between Mn spins, by mapping the spectra of 
magnetic excitations (spin flips) onto a classical Heisen- 
berg Hamiltonian with coupling constants J„, ranging 
from first (n = 1) to sixth (n = 6) nearest neighbors. 
The effects of clustering on the nearest-neighbor pair- 
exchange coupling, J\i have been investigated by examin- 
ing specific (i.e., non-random) configurations with three 
and four Mn atoms in the 128-site supercell: We have 
established that clustering tends to weaken the magni- 
tude of the nearest-neighbor exchange coupling. On the 
other hand, we have found that clustered structures of 
Mn atoms have the lowest total energies, a result which 
may be of importance in a realistic discussion of anneal- 
ing and/or diffusion effects. From calculations on random 
configurations we have also been able to determine the 
behavior of J n with x, for fixed n: in most cases the ex- 
change couplings get weaker as the concentration of Mn 
atoms is increased. This is consistent with the experi- 
mentally observed fact that there is an optimum range 
of Mn concentrations (whose quantitative determination 
requires a careful consideration of other disorder effects) 
in which the critical temperatures are the highest. 

For fixed Mn density, we have found that the calcu- 
lated J n favor a ferromagnetic ground state, and have 
decreasing magnitude as the distance between spins in- 
creases (cf. Table I and Figs. 2 and 4). The non- 
monotonic behavior is attributed to directionality effects; 
by the same token, deviations in the sign of J n were 
found only at large n (=5), when its magnitude is al- 
ready greatly reduced with respect to the nearest neigh- 
bor value. The discrepancy of the present results with 
respect to recent calculations by Xu et al.^ may be 
attributed to the fact that their muffin-tin calculations 
are not full potential; they therefore do not fully repro- 
duce the crucial role played by the directional sp 3 bonds 
and by the hole p-states. Also, due to the quite sig- 
nificant variations of the calculated exchange couplings 
with configurations and Mn concentration, we empha- 
size that estimates of the critical temperature obtained 
via exchange couplings thus obtained are clearly open to 



question. We should also stress that the present results 
corroborate that the Mn-Mn ferromagnetic effective cou- 
pling in Gai-^Mn^As is intermediated by localized holes 
leading to an antiferromagnetic (non-RKKY) coupling of 
each Mn spin, as previously notedj^LS and recently con- 
firmed experimentally. 39 Therefore, the inescapable con- 
clusion is that the main feature of a conventional free- 
electron-like or perturbative RKKY interaction should 
be ruled outs^i^ in the case of Gai_ x Mn x As. 



As a final point, some comments regarding future per- 
spectives are in order. From one side, investigations using 
a similar procedure as employed here (where the disorder 
is explicitly included) of how impurities, such as intersti- 
tial Mn and As anti-sites, alter the effective exchange in- 
teractions are relevant. For a given Mn configuration, it 
should be interesting to see how the results depend on the 
relative position of the defects. On the other hand, our 
results raise some questions whose answers are not com- 
pletely trivial: (i) The fact that the effective exchange 
interactions change with the Mn configuration make it 
clear that the use of a Heisenberg model, at least a simple 
one where only pair-interactions are considered, should 
be viewed with caution. It is not obvious that exten- 
sions of the Heisenberg model to triplets or even larger 
cluster interactions will remedy this fact; (ii) The use of 
ab initio calculations has been very important in order 
to provide a correct picture of the electronic structure 
of these systems. One of its great merits is the possibil- 
ity of obtaining model-free results. However, whenever 
one needs to make predictions about the critical temper- 
ature (T c ), models have to be used. For instance, from 
ab initio results one may extract effective exchange pa- 
rameters, as in the present work, and then via mean field 
or more sophisticated methods, like Monte Carlo calcu- 
lations, it is possible to calculate T c . Two crucial steps in 
this procedure are questionable. The first one is the use 
of a Heisenberg model, as already mentioned. The other 
is the use of a small supercell approximation. Calculat- 
ing the critical temperature via any effective methodol- 
ogy that is based on small supercell ab initio calcula- 
tions, even if this effective approach allows the search 
of a large number of distinct configurations, has a great 
risk of being nonsense, since, as we have shown, the ex- 
change interactions depend sensitively on the Mn distri- 
bution. The root of the above problems is the necessity 
of introducing a model hamiltonian in order to extract 
excited states of the system associated with spin excita- 
tions. A possible solution to this problem could be the 
use of a semi-empirical hamiltonian with a tight-binding 
descrition for the host material coupled with a many- 
body, atomic-like description for the Mn atoms. The 
manyfold of low-energy states representing the different 
Mn spin orientations, that will be obtained upon diago- 
nalization of such a hamiltonian^, will replace the states 
obtained via the effective (but questionable) Heisenberg 
hamiltonian. 
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Appendix 

Here we discuss the case of three Mn atoms, for the 
two disorder configurations displayed in Figs. 1(a) and 
1(b), and chosen according to the SQS algorithm. For 
the configuration in Fig. 1(a) (£ = 4 in Table II), the 
Hamiltonian [see Eq. flfl] may be written, having in mind 
that the interactions with spins on image sites are to be 
taken into account, as 

H = J1S1 • S 2 + J 3 S 2 • S 3 + 4J 6 Si • S 3 (A.l) 

where the absence of second-neighbor, fourth-neighbor, 
and fifth-neighbor interactions should be noticed. 

In a similar way, for the configuration in Fig. 1(b) 
(£ = 9 in Table II), the Hamiltonian is given by 

H = J 3 Si ■ S 2 + 2J 4 Si • S 3 + 2J 5 S 2 • S 3 (A.2) 

where one notes the absence of first-neighbor, second- 
neighbor, and sixth-neighbor interactions. 

One may perform DFT-GGA calculations, and obtain 
the total energies for SQS configurations with all Mn 



5 = 5/2 atoms aligned with each other, as well as for 
an increasing number of flipped Mn total spins. The 
total-energy differences with respect to the aligned states, 
{AE}, may then be obtained via an effective classical 
Heisenberg model with appropriate J n exchange cou- 
plings up to n = 6. 

For the disorder configuration in Fig. 1(a), noticing 
that classically one has Sj • Sj = ±25, it is straightfor- 
ward to obtain, using eq. (A.l), for the total energies 
of configurations with appropriate flipping of Mn total 
spins 

E =f (+Ji + J 3 + 4J 6 ) (A.3) 
El =f (~Ji + J 3 -4J 6 ) (A.4) 
El =f {-Ji- J 3 +4J 6 ) (A.5) 
El =f (+J X - J 3 - 4J 6 ), (A.6) 

where the lower index indicates the number of flipped 
spins (from +5/2 to - 5/2), and the upper index labels 
which spin was flipped. The differences in corresponding 
Heisenberg energies are therefore 

Ai_„ =-f(2J!+8J 6 ) (A.7) 

A2-0 =-f(2J! + 2J 3 ) (A.8) 

A3-0 =_-f (2J 3 +8J 6 ), (A.9) 
and one may thus obtain J\ = —16.2 meV, J 3 = —3.0 
meV, and Je = —1.2 meV from the calculated first prin- 
ciples differences in total energies, i.e., A^q = 264 meV, 
A 2 _ = 241 meV, and A 3 _ = 99 meV. 
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